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SUMMARY 

Applications of CFD methods to determine the regimes of applicability of nonlinear models de- 
scribing the unsteady aerodynamic responses to aircraft flight motions are described The potential 
advantages of computational methods over experimental methods are discussed and the concepts underly- 
ing mathematical modeling are reviewed The economic and conceptual advantages of the modeling pro- 
cedure over coupled, simultaneous solutions of the gasdynamic equations and the vehicle’s kinematic 
equations of motion are discussed The modeling approach, when valid, eliminates the need for costly 
repetitive computation of flowfield solutions For the test cases considered, the aerodynamic model- 
ing approach is shown to be valid 


1 INTRODUCTION 

Predictions of aircraft flight motions, design of aircraft control systems, and development of 
realistic piloted flight simulators all hinge on accurate knowledge of the unsteady aerodynamic 
forces and moments acting on the maneuvering aircraft Prediction of these unsteady airloads is 
complicated by the fact that the instantaneous flowfield surrounding the aircraft 13 not determined 
solely by the instantaneous values of the motion variables such as angle of attack, pitch rate, and 
control deflection angles In general, the instantaneous state of the flowfield, and thus the load- 
ing. depends on the history of the motion, that is, on all states experienced by the flowfield during 
the maneuver prior to the instant in question Time-history effects are accentuated by requirements 
for aircraft agility, which lead to flight in the high-angle-of -attack regime and rapid changes in 
the orientation of the aircraft As a result, highly-maneuverable aircraft can experience nonlinear 
and unsteady airloads due to large regions of three-dimensional separated flow, concentrated vortex 
flows, and the presence and movement of shock waves 

Complete evaluation of time-history effects in wind-tunnel experiments would require the avail- 
ability of an apparatus capable of simulating, at suitable rates, the complete range of motions the 
aircraft would experience in flight, and capable of measuring the aerodynamic response to those mo- 
tions By coupling this apparatus with a computer, the aircraft’s equations of motion could be inte- 
grated in time from specified initial conditions The computed aircraft flight attitude and angular 
rates would determine the new position of the apparatus, while the measured aerodynamic forces and 
moments would be input to update the computed aircraft motion Although, in principle, capable of 
accounting for all time-history effects, such a general wind-tunnel apparatus is, unfortunately, far 
more easily envisioned than constructed As a result the conventional method of predicting aircraft 
flight motions relies on an alternative approach, that of mathematical modeling 


In the modeling approach, one attempts to specify a form for the aerodynamic response that re- 
mains the same in the determination of the aerodynamic response to all motions of interest Ideally, 
a mathematical model enables one to synthesize the response to a general maneuver from the known 
aerodynamic responses of the aircraft to a limited number of specific, characteristic motions Wind- 
tunnel evaluations of the responses to the characteristic motions are, in principle, more easily ac- 
complished than simulations of the general maneuvers The development of mathematical models capable 
of accurately describing the variety of aerodynamic phenomena acting on aircraft maneuvering in the 
nonlinear flight regime has been a subject of ongoing interest (Refs 1,2), and continues to be ad- 
dressed in this symposium (Refs 3,4) 


The utility of a proposed mathematical model is dependent upon the range of vehicle motions and 
aerodynamic phenomena which it can encompass. Determining the range of validity of a candidate math- 
ematical model would require 1) determining the aerodynamic response of an aircraft to the set of 
characteristic motions called for by the model, 2) predicting a wide variety of aircraft motion his- 
tories using the model and the determined responses, and 3) comparing the predicted motion histories 
against actual flight histories The validity of the model would be demonstrated by a match between 
the predicted and flight motion histories While such a model validation procedure is straightfor- 
ward in principle, it is difficult to accomplish in practice Flight-test data are expensive to ob- 
tain, and accurate determination of motion histories during extreme maneuvers is difficult (see Ref 
5 for a discussion of flight-test difficulties) Further, accurate wind-tunnel determination of the 
unsteady responses to the characteristic motions is itself nontrivial (see Refs 6,7 for a descrip- 
tion of relevant experimental techniques) In addition, in the comparison of motion histories, small 
changes in the values of the terms of the model may integrate to produce large deviations between 
the actual and predicted motions Thus, few, if any, satisfactory model validations have been ac- 
complished for a general set of flight maneuvers 

Fortunately, the remarkable advances which have occurred in computer technology and computational 
fluid dynamics now offer the aerodynamic ist a promising alternative approach Today, we can envision 
utilizing CFD techniques to study unsteady three-dimensional flows, both to investigate aerodynamic 
time-history effects and to validate aerodynamic mathematical models 
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A straightforward computational approach to investigating time-history effects (assuming the 
availability of codes and computer resources adequate to solve suitably modeled equations of turbu- 
lent flow around an aircraft) would be to solve the flowfield equations simultaneously with the air- 
craft’s equations of motion Such a coupled approach is shown schematically m Fig la This is. 
of course, the CFD analog of the general wind-tunnel apparatus described above, with computation of 
the instantaneous flowfield replacing the need for measuring the aerodynamic response By using CFD 
techniques one avoids the experimental difficulties which hinder wind-tunnel measurements of unsteady 
flows Results of these coupled computations would be complete time histories of the aerodynamic re- 
sponse and of the vehicle motion It is noteworthy that computations involving the coupled-equations 
approach have recently been carried out for several unsteady two-dimensional inviscid flows (Refs 
8 - 10), and for at least two unsteady two-dimensional viscous flows (Refs 11,12) However, lack 
of computational resources has, to date, precluded undertaking the analogous computations for (even 
steady) three-dimensional viscous flows 

Although coupling of the flowfield equations and the aircraft's inertial equations of motion is, 
in principle, an exact approach to accounting for aerodynamic time-history effects in predicting the 
response to arbitrary maneuvers, it inevitably will be a very costly one This will be particularly 
true for maneuvers at high incidence, where the airloads depend nonlinearly on the motion variables 
Under such conditions the aircraft can experience widely varying motion histories, even if they are 
started from closely-spaced initial conditions Thus, to completely evaluate an aircraft's perfor- 
mance envelope, a large number of computational cases, each involving the coupled equations, would 
be required to cover all possible sets of initial conditions Since the motion and the aerodynamic 
response are inextricably linked in the coupled approach, the flowfields must be recomputed for each 
change in initial conditions 

In contrast to the coupled-equations approach, utilizing mathematical modeling in conjunction 
with CFD techniques can eliminate the need for repetitive computation of flowfield solutions, and 
will lead to more efficient use of computational resources The modeling approach is shown schemati- 
cally m Fig lb As has been previously mentioned, m formulating a model one attempts to identify 
a set of characteristic motions from whose aerodynamic response one can generate the response to gen- 
eral maneuvers Within the regime of validity of the mathematical model, computational evaluation 
of the aerodynamic terms specified by the model would be required only once, whereupon they could be 
utilized over a range of motion variables and flight conditions. Flight motions could then be pre- 
dicted by solving the aircraft's equations of motion independently of the flowfield computations, at 
far less cost 

Presupposing the availability of codes adequate to solve the equations governing turbulent flows 
around an aircraft, evaluation of the validity of a candidate aerodynamic mathematical model can it- 
self be easily accomplished with the use of CFD techniques Two validation procedures can be envi- 
sioned The first method parallels the experimental validation procedure discussed above, that is, 
by comparing vehicle motion histories To carry out the validation, the aerodynamic responses to 
the characteristic motions would be evaluated from flowfield computations Once obtained, these re- 
sponses can be used, together with the model, to predict a series of motion histories On the other 
hand, by using the coupled-equations approach, motion histories which take exact account of time- 
history effects can be obtained for the same set of initial conditions Both series of motion his- 
tones are subject to precisely the same set of approximations made in computing the flowfield Dif- 
ferences, if any, between the motion histories must be attributable to the assumptions made in the 
modeling process Thus, the validity of the mathematical model would be demonstrated by agreement 
between the motion histories The second computational model validation procedure would involve com- 
puting the aerodynamic responses to the characteristic motions, and using the mathematical model and 
the computed responses to obtain the aerodynamic (rather than the motion) response to a complex, but 
specified, motion The response to the complex motion would also be computed directly, and compared 

.to the one obtained using the mathematical model Here the validity of the model would be demon- 
strated by agreement between the two aerodynamic response histories 

In this paper we review two mathematical model validation efforts effected with CFD techniques 
In the first. Chyu and Schiff (Ref 10) utilized two-dimensional unsteady-flow computations to val- 
idate a nonlinear aerodynamic mathematical model for the case of a freely-deflecting flap hinged to 
a stationary airfoil and immersed in a transonic flow Validation of the model was demonstrated by 
close agreement of flap motion histories obtained with the model and those obtained with the coupled- 
equations approach More recently, Katz and Schiff (Ref 13) investigated the validity of a non- 
linear aerodynamic mathematical model for low-rate multi-degree-of -freedom motions of a delta wing 
maneuvering at large angles of attack. Although computation of three-dimensional unsteady viscous 
flows by means of modeled equations of turbulent flow is not yet feasible, computation of unsteady 
•three -dimensional inviscid flows by means of potential equations is currently possible The authors 
of Ref 13 employed a nonlinear vortex-lattice method to compute the aerodynamic responses to the 
characteristic motions called for by the model Aerodynamic responses to complex specified motions 
were generated by means of the model and also from direct flowfield computations The validity of 
the mathematical modeling concepts for the delta wing maneuvering in the high-angle-of-attack regime 
was demonstrated by close agreement of the corresponding force histories In the following sections 
we discuss the procedures and results of Refs 10 and 13, focusing in particular on the valida- 
tion of the mathematical models and on the costs of the modeling approach relative to those of the 
coupled- equations approach 


2. AERODYNAMIC MATHEMATICAL MODELS 

The form that an aerodynamic mathematical model takes depends upon the coordinates used and vari- 
ables used to describe the motion, and on the types of aerodynamic phenomena that the model can en- 
compass In a series of papers, Tobak and his colleagues used the concept of a nonlinear indicial 
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function to derive aerodynamic models to describe, within certain assumptions, nonplanar maneuvers 
of bodies of revolution (Ref 14) and of aircraft (Refs 15,16) In this section an alternative 
derivation is reviewed which puts clearly into perspective the assumptions made in deriving the mod- 
els, and the resulting limitations of the resulting models to accommodate certain types of aerody- 
namic phenomena 

2 1 Model Describing Aircraft Maneuvers 

Consider a maneuvering aircraft whose center of mass travels an essentially straight-line path at 
constant flight velocity The maneuvers can be described in terms of an aerodynamic axis system, in 
which the orientation of the aircraft relative to the oncoming wind is specified in terms of pitch 
and roll coordinates (Fig 2) In this axis system the resultant angle of attack, <r, is defined as 
the angle between the body-fixed longitudinal X axis and the flight velocity vector The plane con- 
taining a is called the reoultant-angle-of -attack plane The roll, or bank, angle 0 is the angle be- 
tween the normal to the resultant-angle-of -attack plane and a body-fixed axis (normal to the X axis) 
which lies in the plane of the wing The resultant-angle-of -attack plane is free to rotate about the 
velocity vector with a rotation rate, measured relative to an inertia frame, of 0 Knowledge of the 
history of the maneuver is equivalent to knowledge of the histories of <r, 0, and 0 

The aircraft undergoes a maneuver beginning at time zero The instantaneous aerodynamic response 
to the maneuver at time t depends, in general, on the entire history of the maneuver That is, the 
response depends on the histories of the motion variables over the interval from 0 to t Mathemati- 
cally, the response is a functional of the motion variables Thus the pitching-moment coefficient, for 
example, can be expressed in the square-bracket notation of Volterra, as 

c m (t) = c m |ff(fl,0(O.*(€)l (i) 


where f is a dummy time variable ranging from 0 to t Note that 0(f) appears in Eq (1) , rather 
than 0(f). because the aerodynamic response is independent of the orientation of the angle-of -attack 
plane (m contrast to that of the pilot) and depends only on the rate at which the o plane is rotat- 
ing 


If the motion is assumed to be differentiable, one can expand each of the variables o(f), 0(f), 
and 0(f) m Taylor series about f = t Knowledge of the past history of the motion is equivalent to 
knowing all of the coefficients of the series Substituting these coefficients for the motion his- 
tories in Eq (1) converts the pitching-moment response from a functional to a function, albeit a 
function of an infinite number of variables Thus, the pitching-moment coefficient becomes 

C m (t) = C m (o{t),o(t),o(t), ,<j>(t),4>(t), ) (2) 

The assumption that the motion history is differentiable immediately restricts the model to cases 
where the aerodynamic forces and moments are continuous If sudden jumps in the aerodynamic re- 
sponse occurred, the aircraft would experience discontinuous accelerations, and this would violate 
the assumption of a differentiable motion history Restriction to continuous forces rules out the 
capacity of the resulting mathematical model to treat the bi-valued aerodynamic responses charac- 
terizing static hysteresis and both time-invariant and time -dependent subcritical bifurcations (cf 
Ref 3) Within the restriction to continuous, single-valued, aerodynamic responses, one can make 
a further assumption that, for the slowly-varying notions typifying aircraft rigid-body motions, the 
aerodynamic response will have only a negligible dependence on the higher-order rate terms such as 
<?(/), cr(f ) , 0(£), 0(£), 0(<) , etc , and thus they may be neglected in Eq (2) This reduces the pitching- 
moment response to a function of a finite, rather than an infinite, number of variables Consistent 
with the assumption of slow motion, one is justified in expanding Eq (2) in a Taylor series about 
<7 = 0, 0 = 0, and 0 = 0 The resulting model is 


C m (t) =C m (o(t),V>(0>0, 0,0) + 

+ O.o.o.o) + — c mi (<r(e) o, o, o) 


(3) 


where 6 is the wing span, and the zeros refer to the motion rates Equivalent expressions for the 
yawing- and rolling-moment coefficients, C n and Ci , and for the axial-, side-, and normal-force co- 
efficients, Cx . Cy . C/v, are obtained by substituting these coefficients wherever C m appears in 
Eq (3) Again, the model applies to slowly varying motions of the aircraft, although the val- 
ues of o and 0 may be large The mathematical model is seen to contain four terms, and each term 
can be identified with a specific, characteristic notion from which it may be evaluated Thus, 
C m [a(t) n 0(0* 0,0,0) is the pitching-moment coefficient that would be evaluated in a steady flow with o 
and 0 held fixed at a(£), 0(0 The second term, Cm 4 . is the contribution to the pitching-moment co- 
efficient due to rolling motion, and can be evaluated for small-amplitude oscillations in 0 about 0 
= const, with o held fixed and 0 fixed at zero Similarly, the term C m<f is the contribution to the 
pitching-moment coefficient due to pitching motions, and can be evaluated for small-amplitude pla- 
nar oscillations in o about a = const, with 0 held fixed and 0 fixed at zero. The last term, C m ^ , 
is the rate of change of the pitching-moment coefficient with coning-rate parameter, 06/21/ , evalu- 
ated at 0 ■ 0, that would be determined from a steady coning motion with o - const, 0 * const, and 0 
= const These characteristic motions are illustrated in Fig 2 
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2 2 Model Describing Flap Motions 

Derivation of an aerodynamic mathematical model applicable to slowly-varying motions of a flap 
on an airfoil parallels that of the model. Eq (3). describing the response to aircraft maneuvers 

For the airfoil flap motions considered in Ref 10, the coordinates illustrated in Fig 3 were em- 

ployed The airfoil was a NACA 64A010 section hinged at the 75% chord point The forward portion of 
the airfoil was held fixed parallel to the oncoming stream, while the flap was free to move about its 
pivot point Positive values of the flap deflection angle, o/ t and of the flap hinge-moment coeffi- 
cient, Ch . are shown in Fig 3 As illustrated, a positive value of Ch would tend to increase the 
flap deflection angle For the single-degree-of -freedom flap motion, the hinge-moment coefficient 
is, in general, a functional of the flap deflection angle history, i a . C*(f) = C*[tf/(£)] By means 
of the same logic described above for aircraft maneuvers, the functional can be reduced to the form 

c h (t) = C h (<T,(t),0) + °-£c h ' t {o f {t), 0) (4) 


where the reference length / is the chord length. As before, the model applies to slowly-varying mo- 
tions of the flap, although the values of the flap deflection angle can be large Again, each of the 
terms in Eq (4) is identified with a particular characteristic motion from which it may be evalu- 
ated Thus, the term 0) is the hinge-moment coefficient that would be evaluated in a steady 

flow with the flap deflection angle held fixed at o/(t) The remaining term, Ck„ f . is the contribu- 
tion to the hinge-moment coefficient due to flap motions, and can be evaluated for small-amplitude 
oscillations of the flap about a mean value of cry held fixed at the instantaneous value of cry(f) 

2 3 Issues of Model Applicability 

Consistent with the assumptions made m their development, aerodynamic mathematical models at the 
level of Eqs (3) and (4) are subject to the following restrictions 1) the response to a steady 
motion is itself steady, 2) the response is a single-valued (although allowably nonlinear) function 
of the orientation of the body, and 3) the responses are linear in the motion rates Restriction 1 
rules out the possibility of modeling time-dependent aerodynamic bifurcations, that is, development 
of time-varying (periodic, quasi-penodic , or chaotic) responses to a steady motion, e g , the time- 
varying flow observed surrounding a stationary circular cylinder in crossflow at Re > 50 However, 
recent modeling efforts (see Ref 3) , based on the concept of FrSchet differentiability of the aero- 
dynamic response, include treatment of such time-varying flows Restriction 2 precludes modeling the 
bi-valued aerodynamic responses characterizing static aerodynamic hysteresis which have been observed 
in cases of vortex asymmetry on slender bodies of revolution and on slender delta wings However, 
this behavior can also be accomodated (see Ref 1. and more recently. Ref 3) The third restric- 
tion precludes modeling the nonlinear dependence of the aerodynamic response on the motion rates, and 
in particular, nonlinear dependence on the coning rate Such nonlinear variation with 4> has been 
observed experimentally (see Ref 17) , and can be incorporated in the model by performing the Taylor 

series expansion about o = 0. tJj = 0, and <j> — <t>[t) rather than about o = rj) = ^ = 0, as was done to obtain 
Eq (3) This would result in an aerodynamic model equivalent to the one derived previously (Ref 
16) using the nonlinear indicial function approach 

The utility of the aerodynamic modeling approach depends on the ability of the model to treat the 
aerodynamic phenomena which occur in flight In applying models analogous to Eq (3) , the general 
aircraft motion is decomposed into a sum of characteristic motions The aerodynamic response to the 
general motion is modeled as a sum of responses to the characteristic motions The actual response 
to the general motion will differ from the modeled response if aerodynamic phenomena excluded in de- 
veloping the model are present The assumption made in developing Eq (3), that the aerodynamic re- 
sponses are continuous, single-valued functions of the motion variables, restricts the model to cases 
where neither hysteresis nor time -dependent aerodynamic bifurcations occur Within this restriction, 
the remaining causes for failure of the model to predict a general response would be either 1) sig- 
nificant nonlinear dependence of the aerodynamic responses on rates of motion within the range of 
rates actually experienced in flight, or 2) presence of significant interactions between responses 
to pairs characteristic motions Examples of such interactions include those between responses to 
pitch oscillations and coning motion or between responses to roll oscillations and coning motion, 

1 e , terms such as or which have been excluded in deriving Eq (3) 

2 4 CFD Requirements for Validating Models 

As discussed in Sec 1, validation of a candidate mathematical model with the use of CFD methods 
involves computation of both the aerodynamic response to characteristic motions and the response to 
either specified complex motions or coupled-equations responses From a CFD standpoint, the boundary 
conditions for the low-rate characteristic pitching, rolling, and coning motions are linear perturba- 
tions in < 7 , 0, and <f> t respectively, about the zero values associated with the steady characteristic 
motion If flowfields were governed by linear partial differential equations, the linear boundary 
conditions guarantee that the responses to the characteristic motions would be linear in the rates 
Also, the response to any general motion could be obtained from superposition of the responses to the 
characteristic motions, for such a solution would satisfy both the PDE's and BC's exactly In this 
circumstance, models at the level of Eq (3) would be exact for all cases 

However, the actual equations governing fluid flow are nonlinear, and the presence of nonlinear 
convection terms in the momentum equations give rise to the possibility of all the aerodynamic phe- 
nomena discussed above Thus, to be useful in assessing the validity of a model, the computational 
method must be based on nonlinear flowfield equations Approximating the flowfield equations with 
linearized PDE's will clearly be unsatisfactory for validating nonlinear models Such models will 
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always appear to be valid, since the linearized flowfield equations cannot capture nonlinear aerody- 
namic effects Thus , validation of limited aerodynamic models analogous to those of Eqs (3) and 
(4) would require, at a minimum, a CFD method capable of demonstrating nonlinear aerodynamic behavior 
with the motion rates In Ref 10. Chyu and Schiff used a time-accurate method based on the non- 
linear Euler equations to demonstrate the validity of the model, Eq (4). for the case of an oscil- 
lating flap The Euler equations are capable of capturing nonlinear behavior of both the steady and 
unsteady flowfield responses Similarly, in Ref 13, Katz and Schiff utilized a nonlinear vortex- 
lattice method to investigate the ability of the equivalent model, Eq (3). to describe the nonlin- 
ear aerodynamic response to low-rate maneuvers of a delta wing In the vortex-lattice technique, the 
flowfield is assumed to be governed by the incompressible potential equation, V 2 <£ = 0 Leading-edge 
separation is modeled by permitting vortex panels to shed from the leading edge of the wing, and the 
free vortex sheet is allowed to roll up in response to the local velocity fields Although the PDE 
which describes the flowfield is linear, the boundary conditions are not Allowing the free vortex 
sheets to roll up permits the method to encompass a nonlinear variation, if present, of the aerody- 
namic response with the rates (see Ref 18) 

Methods based on nonlinear, inviscid flowfield equations would appear adequate to validate aero- 
dynamic models analogous to those of Eqs (3) and (4) These computational methods cannot, however, 
incorporate aerodynamic phenomena such as onset of vortex asymmetry, vortex breakdown, and time- 
varying response to fixed boundary conditions Thus, validation of aerodynamic models which purport 
to include these phenomena will require the use of methods based on nonlinear viscous flowfield equa- 
tions 


3 VALIDATION FOR MANEUVERING DELTA WING 

In the study reported in Ref 13, CFD solutions were carried out for the unsteady three- 
dimensional flowfield surrounding a sharp-leading-edge delta wing to demonstrate the validity of the 
multi-degree-of -freedom mathematical model, Eq (3) The flowfield was assumed to be governed by 
the potential equations, and a vortex- lattice method (VLM) was applied to solve the time-dependent 
equations 

3 1 Model Validation Procedure 

The particular wing considered in Ref 13 was a slender delta wing having an aspect ratio of 
unity (leading-edge sweep angle = 7596°) The center of mass was fixed at the wing half-chord (x/c 
= 050) The values of the resultant angle of attack and roll angle considered are shown in Fig 4 
The resultant angle of attack a ranged from 20° to 30°. while the roll angle ^ ranged from 0° (wing 
level) to 10° The dimensionless rates ob/2U . tpb/2U , and <£6/2(7. ranged up to 0 15 The procedure 
utilized for validating the aerodynamic mathematical model applicable td the wing had three phases 
These were 

1) Evaluate the aerodynamic responses to the characteristic motions from vortex-lattice computa- 
tions 

2) Generate aerodynamic force and moment response histories to prescribed complex motions using 
the aerodynamic mathematical model and the aerodynamic data evaluated in phase 1 

3) Compare the histories obtained in phase 2 with force and moment histories that are, in princi- 
ple, exact within the framework of the computational technique, namely those obtained by directly the 
vortex-lattice method to compute the response to the identical complex motions 

Demonstration of the validity of the mathematical modeling approach, as applied to the delta wing 
maneuvering in the high-angle-of attack regime, hinged on finding close agreement between the force 
and moment time-histories obtained from the two approaches 

3 2 Numerical Technique 

In the computational procedure, the surface of the wing is divided into a number of bound vortex 
panels The strengths of the bound panels are determined at each time step during the computation 
to enforce the boundary condition that there be no flow through the solid wing Information describ- 
ing the wing’s maneuver enters the computation through the solid-surface boundary condition Time- 
evolution of the wake behind the wing is modeled by allowing vortex panels to shed from the trailing 
edge at each time step These wake panels have fixed strength and, upon leaving the wing, move with 
the local fluid velocity 

When a delta wing is maneuvering at high angle of attack, flow separates near the wing leading 
edges, and the separated fluid rolls up above and behind the wing to form concentrated vortices For 
a sharp-edged wing, the separation line is essentially fixed at the sharp leading edge, and does 
not vary with changes in Reynolds number In the vortex-lattice method, leading-edge separation 
is modeled in a manner analogous to that of trailing-edge separation, by allowing vortex panels to 
shed from lines of separation specified to remain at the leading edges and permitting them to move 
with the local flow velocity Specification of the bound vortex strengths and free vortex posi- 
tions yields the pressure distribution on the wing and, in turn, the nonlinear unsteady airloads (cf. 
Refs 18-21). Details of the numerical method and a discussion of the accuracy of the computed re- 
sults are found in Ref. 18 
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3 3 Aerodynamic Response to Characteristic Motions 

To obtain the aerodynamic data required by the mathematical model. Eq (3), computations were 
carried out for the wing in each of the four characteristic motions shown in Fig 2 at each value of 
resultant angle of attack and roll angle shown in Fig 4 The steady-state term Cfc(a(t)i 0(0) » where 
Ci t denotes any of the force or moment coefficients, was obtained from a computation in which the re- 
sultant angle of attack and roll angle were held fixed, and the flowfield was allowed to evolve un- 
til it reached a steady state In an analogous manner the term 0(f)) was obtained from a se- 

ries of computations for steady coning motion in which the resultant angle of attack, roll angle, and 
coning-rate parameter were fixed, and the flowfield was allowed to evolve to a steady state Note 
that to an observer fixed in the moving wing, the flowfield due to a steady coning motion is indeed 
time-invariant The coefficient was then determined from the observed rate of change of the moment 

with coning-rate parameter. dC^ld[4>bl2U ) , evaluated at <f> * 0 

The aerodynamic coefficient due to pitch oscillations, Cfc»(<r(t)i0(O) . was evaluated from small- 
amplitude harmonic pitch oscillations about the mean values of resultant angle of attack and roll an- 
gle shown in Fig 4 The wing was specified to move according to 

a = a 0 4- o\ smcjj t 

0 = 0o (5) 

<t> = 0 


The amplitude of the harmonic motion, o \ , was specified to be less than 2° The aerodynamic damping 
coefficient was evaluated from the component of the aerodynamic response that was 90° out of phase 
with the wing motion This rationale is easily seen by substituting the conditions describing the 
pitch oscillations (Eq (5)) into the aerodynamic model to obtain (after a Taylor-series expansion 
about a ~ cq and omission of terms of 0(aj)) 


c*(i) = Cfc(po,t/i 0 ) + ffi V °) 4- cosull tC ki (a 0 ,tl> 0 ) 


(6) 


The coefficients in Eq (6) were obtained from a Fourier integration of the response over one cycle 
of the motion, as shown for the normal-force coefficient response in Fig 5. Thus 


i r 2K 

Cfc(a o ,0 o) = ^ J o c k(t)d(uit) 

(7) 

i r 2 * 

Ck„{o o,0o) = / Cfc(t) sinwifrf(u»i<) 

(8) 

Jo 

2U f 2 * 

CkA°o,il>o) - lrbu>l( , l J <?*(«) coscditd(wit) 

(9) 


The steady-state coefficient and its slope can be obtained either from the computations of the os- 
cillatory motion (Eqs (7) and (8)) or, preferably, from the computations of the steady motion de- 
scribed earlier 


In an analogous manner, the coefficient due to roll oscillations. Cjt^ (cr (f ) , t/>(*)) . was evaluated for 
small-amplitude harmonic roll motions where 


a — Co 

0 = 0o 4- 0i smu> 2 < (10) 

0 = 0 


and the amplitude of the motion, 0i , was specified to be less than 2° The roll damping coefficient 
was obtained from 


Ck+(c r o ,0o) 


2U 

7t6cJ2 0I 



Ck(t) cos U 2 td(cj 2 f) 


( 11 ) 


The results of the computations for the characteristic motion are shown in Figs 6-8 Generation of 
these diagrams of the aerodynamic coefficients required 36 individual computations, one for each of 
the four characteristic motions at the nine combinations of resultant angle of attack and roll angle 
shown in Fig 4 The results confirm that, over the range of angles of attack and roll and over the 
range of rates considered, the aerodynamic responses to the characteristic motions are linear in the 
motion rates, and are single-valued functions of the angles Thus, barring the presence of signifi- 
cant nonlinear interactions between the responses to pairs of characteristic motions, the aerodynamic 
model, Eq (3), should prove to be valid 
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3 4 Aerodynamic Response to Prescribed Complex Motions 

The prescribed complex motions combined pitch oscillations, roll oscillations, and coning motion 
The combined motions all had the basic form 


a ~ Co + o\ sin a - w\ 0 \ cos uM 

0 = V >0 + Ipl sinu> 2 <, 0 = u /2 01 cosu >2 1 ( 12 ) 

<i> = <t> - W 3 

Aerodynamic response histories of the pitching-moment, rolling-moment, and normal-force coefficients 
were computed from Eq (3), with the aerodynamic coefficients obtained from table look-ups of the 
data shown in Figs 6-8, and values of <7, 0, and <j> obtained from Eq (12) The aerodynamic re- 

sponses to the combined motions were also obtained from direct VLM computations As was mentioned 
earlier, the use of the identical vortex-lattice method to evaluate both the nonlinear responses to 
the characteristic motions and the responses to the general large -amplitude motions ensures a consis- 
tent treatment of the time-history effects Thus discrepancies, if present, between results obtained 
using the modeling approach and those obtained from the direct VLM computations must be attributed to 
the inadequacy of the aerodynamic model 

3 5 Validation of the Mathematical Model i 

The motions considered in Figs 9 and 10 combine pitch and roll oscillations, having ampli- 
tudes of 3° and reduced frequencies ranging up to 0 15. superimposed on a steady coning motion in 
which o 0 = 25° and 0 O = 5° Typical responses of the pitching-moment, rolling-moment, and normal- 
force coefficients to these motions are shown in Figs 9 and 10 as functions of wing-chord lengths 
of travel. Ut/c In each figure the dotted lines show the histories obtained from the mathematical 
model, Eq (3), while the solid lines indicate the results obtained from direct VLM computations 
For the direct computations, the overshoot indicated at the beginning of each time-history occurs be- 
cause the motion was started impulsively from rest It will be recalled that the mathematical mod- 

els described m Sec 2 were obtained under the assumption of slowly-varying motions, and are not 
expected to model the impulsive start If the short initial transient period is excluded, in both 
cases the aerodynamic responses obtained from the model show reasonable agreement with those obtained 
from direct computation The differences between the results obtained from the model and those from 
direct VLM computation can be attributed to errors m the interpolations in the table look-ups based 
on Figs 6-8 that were required to obtain the coefficients of the model This would indicate that, 
for these cases, no significant interactions existed between responses to pairs of characteristic mo- 
tions. such as the interaction between the responses to the pitch oscillation and coning motion, or 
between the responses to the roll oscillation and coning motions 

The agreement shown in Figs 9 and 10 would tend to indicate that, over the range of angles and 
motion rates considered in Ref 13, the mathematical model presented in Eq (3) is adequate to de- 
scribe the aerodynamic response to complex motions of the delta wing In actuality, for the range of 
pitch and roll rates considered, the contributions of the pitch-damping and roll-damping terms in Eq 
(3) were almost negligible As a result, the aerodynamic interactions between pairs of the charac- 
teristic motions, which are of higher order than the damping terms, must be negligible In these 
circumstances it is not surprising that the mathematical model would appear to be validated Thus, 
while the VLM computations confirm the validity of the model for the cases considered, the cases 
themselves do not conclusively demonstrate the limits of the range of motions and rates for which the 
model is valid Nevertheless, the procedure discussed indicates the way in which CFD methods can be 
used to validate a candidate mathematical model 


4 VALIDATION FOR FLAP MOTIONS 

In the study reported in Ref 10, CFD solutions were obtained for the two-dimensional unsteady 
transonic flowfield surrounding a wing and moving flap (Fig 3) to demonstrate the validity of the 

nonlinear mathematical modeling concepts The flowfield was assumed to be governed by the time- 

dependent mviscid Euler equations, and an implicit time-accurate finite-difference technique was ap- 
plied to obtain the solutions The flowfield equations were solved on a moving body-conforming com- 
putational mesh, which deformed in response to the flap motion (Fig 11) Details of the numerical 

technique, boundary conditions, and mesh-generation method are found m Ref 10 

4 1 Model Validation Procedure 

The airfoil considered in Ref 10 was a NACA 64A010 section, hinged at the 75% chord point The 
forward portion of the airfoil was held fixed parallel to the oncoming stream, while the flap was 
free to pivot about its hinge point with deflection angles ranging up to 20° The free-stream Mach 
number was held fixed at 0 8 Validation of the aerodynamic mathematical model involved 

1) evaluating the aerodynamic response to the characteristic motions from Euler-equation computa- 
tions 

2) utilizing the equation governing mechanically unconstrained motions of the flap, 

Ib f {t)=qSlC h (t) (13) 


together with the aerodynamic mathematical model, Eq (4). to predict a series of flap motion histo- 
nes 
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3) simultaneously solving the flowfield equations and flap-motion equation to obtain "exact" flap 
motion histories with the same initial conditions 

The validity of the modeling approach for this nonlinear, transonic flow case was demonstrated by 
close agreement of the motion histories obtained from the two approaches 

4 2 Hinge-Moment Response to Characteristic Motions 

To obtain the hinge-moment coefficients required by the mathematical model, computations were 
carried out for the flap in both required characteristic motions The steady-state term C^(of) was 
obtained from a computation in which the flap deflection angle was held fixed A time-invariant mesh 
was employed, and the flowfield was allowed to evolve until it reached a steady state The resulting 
static hinge-moment coefficient, obtained from spatial integration of the surface-pressure distribu- 
tions, is shown in Fig 12 The hinge moment is statically stabilizing, tending to oppose the flap 
deflection The nonlinear behavior of the static hinge moment is associated with rearward movement 
and Increase in strength of the upper-surface shock wave with increasing flap deflection For low 
flap deflection angles, at the Mach number considered, the shock wave was located ahead of the flap 
hinge point As the flap deflection was increased beyond 3°, the shock moved onto the flap, and its 
rearward movement along the flap, for 3° < Of < 10°, caused the greatly increased slope of the hinge- 
moment curve For values of Of > 17°, the shock wave was essentially fixed at the trailing edge 

The contribution to the hinge-moment coefficient due to flap oscillations, Ch n/ . were evaluated 
from the periodic response to small-amplitude harmonic flap oscillations The flap motion was speci- 
fied as 


<V(0 = °U +°h smwt 


(14) 


where Of ia ranged from 0° to 20°, and oj x was small, usually 0 5°. The reduced frequency of the mo- 
tion. k = ul/V , was held fixed at 0 15 for all cases Starting from an initial steady solution ob- 
tained at a/ computations were carried out for three cycles of the motion to ensure that a pe- 

riodic solution had been obtained The hinge-moment damping coefficient was evaluated from the com- 
ponent of the aerodynamic response that was 90° out of phase with the flap motion The results are 
shown in Fig 13 as a function of the mean flap deflection angle Note that at the transonic flow 
conditions considered in Ref 10, the damping coefficient is a highly nonlinear function of the flap 
deflection For values of <7/ m < 3° the coefficient is negative, or dynamically stabilizing However, 
for mean flap deflections ranging between 3° and 17°, the coefficient is positive (dynamically desta- 
bilizing) and would cause an unconstrained flap oscillation to increase in amplitude The nonlin- 
ear behavior of the damping coefficient is qualitatively related to the behavior of the static hinge 
moment In particular, the decrease in dynamic stability is linked to the increase m slope of the 
static hinge-moment curve (see Ref 22 for a detailed discussion) 

4 3 Flap Motion Histories 

Oscillatory time-histories of the flap motion were generated using the flap equation of motion, 

Eq (13), with the instantaneous hinge-moment coefficient specified by the nonlinear aerodynamic 
model In these computations, the flap moment of inertia was chosen to give a value of reduced fre- 
quency close to the one specified for the characteristic-motion computations After initial values 
of the flap deflection angle and velocity were specified, the equation of motion was solved numer- 
ically to obtain the motion histories At each time step the hinge-moment coefficient was speci- 
fied by Eq (4), where the terms and (*/(*)) were obtained from table look-ups in Figs 

12 and 13, respectively Corresponding histories of the flap motion were also predicted using the 
coupled-equations approach 

Time-histories of flap motion, generated with both the modeling approach and the coupled- 
equations approach are shown in Figs 14 and 15 The motion resulting when the flap was released 
from rest with an initial deflection angle o/ Q =40° is shown in Fig 14 In this case the stabi- 
lizing portion of the hinge-moment damping curve governed the motion, and the amplitude of the os- 
cillation decayed smoothly In contrast, when the flap was released with a slightly larger initial 
deflection, Of n =45°, the dynamically destabilizing portion of the damping curve caused the amplitude 
of the oscillation to grow rapidly (Fig 15) In both cases, however, the motion histories obtained 
using the aerodynamic mathematical model were in good agreement with the "exact" motion histories ob- 
tained from the coupled equations. This confirmed the ability of the mathematical model, Eq (4), 
to describe the unsteady aerodynamic response in this highly nonlinear transonic flow On the other 
hand, a motion history generated with an aerodynamic model which does not account for nonlinear un- 
steady aerodynamic contributions (i e., Eq (4), but with the value of the damping term held 

fixed for all o/ at the value obtained at oj =0°) failed to predict the undamped growth of the of the 
flap oscillation (Fig 16) 


6 DISCUSSION 

The agreement between the force and moment histories obtained for the delta wing using the mod- 
eling approach and those obtained from direct CFD computations confirmed the ability of the model to 
describe the unsteady aerodynamic response to complex motions, at least within the limited range of 
motions and rates considered Similarly, the agreement between the flap motion histones obtained 
with the coupled-equations approach and those generated with the model confirmed the capacity of the 
model to adequately describe the unsteady, nonlinear aerodynamic response to arbitrary flap motions 
The success of the modeling approach in this case points the way to approach future problems involv- 
ing unsteady motions In our view, for cases where mathematical models are valid, the modeling ap- 
proach will be more economical and more versatile than the coupled-equations approach Further, use 
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of the modeling approach will give better insight into the underlying physics than can be obtained 
through use of the coupled-equations approach 

First, regarding costs, we note that the computation of motion histones from the vehicle equa- 
tions of motion requires negligible computational effort in comparison with that needed for computa- 
tions of time-dependent flowfields Consequently, with the modeling approach, once the initial ef- 
fort of evaluating the aerodynamic responses to the characteristic motions is expended, computation 
of motion histones would be relatively inexpensive Regarding versatility, the modeling approach 
makes it easy to introduce changes into the aircraft's equations of motion (e g , changes in vehicle 
mass or moment of inertia, or inclusion of a model of a control system) and to evaluate their effects 
at low cost, since the aerodynamic data within the mathematical model will remain unchanged In con- 
trast, in the coupled-equations approach the simplest change in the aircraft's equations of motion 
would require a complete reevaluation of the flowfield and motion response 


Second, the modeling approach would appear to give better insight into the physics governing the 
unsteady flow than would the coupled- equations approach If an undamped or divergent motion results 
from coupled-equations computations, it would be difficult to identify the aerodynamic phenomena 
causing the instability On the other hand, computations carried out in terms of the characteris- 
tic motions permit an investigation of the underlying aerodynamic mechanisms For the transonic flap 
motions considered in Ref 10, these computations indicated that it was the rearward movement of the 
upper-surface shock wave that caused the large change on the slope of the static hinge-moment curve 
It was also possible to show how the change in the static hinge-moment coefficient was related to 
the destabilizing behavior of the hinge-moment damping coefficient Further, the modeling approach 
is compatible with established methods for determining the stability of motions In the case of the 
flap, knowledge of the behavior of the damping coefficient with increasing deflection, obtained from 
the characteristic motion computations, permitted prediction of the change from a damped to an un- 
damped oscillation with increased initial deflection angle that was subsequently observed Thus, the 
nonlinear modeling approach would appear to be the method of choice in the design of flight control 
systems and m flight simulations 


6 CONCLUDING REMARKS 


Applications of CFD methods to demonstrate the regimes of validity of nonlinear aerodynamic math- 
ematical models for the case of a maneuvering delta wing, and for the case of an oscillating flap on 
an airfoil in transonic flow, are discussed The assumptions underlying the development of the mod- 
els are reviewed, and the ensuing limitations on the aerodynamic phenomena that the resulting mod- 
els can accommodate are discussed The class of mathematical models considered, where the nonlinear 
aerodynamic responses are continuous single-valued functions of the motion variables, was found to 
be adequate for the types of aerodynamic phenomena that were in play Extension of the mathematical 
models to encompass a wider variety of possible aerodynamic phenomena, including bi-valued aerody- 
namic responses and time-varying responses to steady motions is currently being pursued 


The economic and conceptual advantages of the modeling approach over that of coupled, simultane- 
ous, solutions of the flowfield and kinematic equations of motion in predicting flight vehicle motion 
histories are illustrated The modeling approach, when valid, eliminates the need for costly repeti- 
tive computation of flowfield solutions when multiple motion histories must be determined 
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Fig. 1. Approaches to applying CFD methods to predict aerodynamic responses and aircraft motion 

histories. 
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Fig 2. Aerodynamic axis system and characteristic motions obtained assuming linear variation of the 

response on the motion rates 
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Fig. 3 Coordinates and notation for flap motion. 



Fig. 4 Range of angle of attack a and roll angle ip considered in determining aerodynamic coeffi- 
cients for the mathematical model (Ref. 13) 
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Fig. 11. Typical body-conforming grids with flap at extreme deflection angles (Ref. 10). 
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